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Condensate-profile asymmetry of a boson mixture in a disk-shaped harmonic trap 
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A mixture of two types of hard-sphere bosons in a disk-shaped harmonic trap is studied through 
path-integral quantum Monte Carlo simulation at low temperature. We find that the system can 
undergo a phase transition to break the spatial symmetry of the model Hamiltonian when some 
of the model parameters are varied. The nature of such a phase transition is analyzed through 
the particle distributions and angular correlation functions. Comparisons are made between our 
calculations and the available mean-field results on similar models. Possible future experiments are 
suggested to verify our findings. 

PACS numbers: 03.75.Mn, 05.30.Jp, 02.70.Ss 



Mixtures of different cold atoms are studied inten- 
sively by both theorists and experimentalists after the 
first experimental observation of a stable double conden- 
sate formed from two different spin states of 87 Rb [l| . In- 
teresting phenomena, such as the formation of molecules, 
can be directly realized in these mixtures 0, Q ■ Further- 
more these mixtures form a rich set of systems that can 
test many-body theory and novel ideas at extremely low 
temperature and quantum limits 0, 0, ■ 

One of the intriguing features predicted to occur in 
such systems, based on the solutions of the Gross- 
Pitaevskii equation 0, H H H EJ H 13 Q, is the 
broken spatial symmetry of the Hamiltonian in the con- 
densate profiles of individual species in the system. This 
is purely a quantum phenomenon because it is not ex- 
pected in the particle distributions of the corresponding 
classical systems. It is known that a double condensate 
can form an inner core from one species and an outer shell 
from the other in a spherically symmetric trap (lH , ll(i| , 
which is a result of the competition between correlation 
and coherence presented in the system [l7j , while preserv- 
ing the symmetry of the Hamiltonian. But this structure 
may become unstable for certain interaction strength and 
total number of particles and the system can undergo a 
quantum phase transition to a state with a broken spa- 
tial symmetry 0,0 — a double-condensate mixture in an 
anisotropic trap can develop an asymmetry in condensate 
profiles along the direction of the weakest confinement, 
according to the mean-field studies H H [H El E . 
But it is not entirely certain whether these predications 
are genuine results of the model system or some artifacts 
of the mean-field approximation. 

In this Letter we report key results from our recent 
study of the two-component hard-sphere boson mixtures 
in a disk-shaped harmonic trap through path-integral 
quantum Monte Carlo simulation, which is exact within a 
controllable sampling variance. This method allows us to 
study many-body quantum systems at finite temperature 
accurately and to make precise predications about the 
model systems studied. The technique has been applied 
successfully to boson cluster in various traps 18] and to 
the two-component system in a harmonic trap that is 



spherically symmetric |l7j . For a double condensate in a 
disk-shaped harmonic trap, we find that each component 
can undergo a phase transition and break the spatial sym- 
metry of the Hamiltonian (or that of the trap potential) 
when certain parameters of the system, such as the inter- 
action range, or the aspect ratio of the axial confinement 
along the z direction and the radial confinement in the 
xy plane, are varied. The results are in qualitativ e ag ree- 
ment with the mean-field studies 0, 0, Q, 0, H3 | 
of this model. In order to understand this symmetry- 
breaking phenomenon in detail, we analyze the influences 
of the interaction range, temperature, and number of par- 
ticles of each species, external potential and elucidate the 
nature of this phase transition. 

Here we consider a two-component mixture of hard- 
sphere bosons described by the model Hamiltonian 
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is the Hamiltonian of the corresponding noninteracting 
system under the trapping potential f7;(r) — miUj\(x 2 + 
y 2 + X 2 z 2 )/2. The parameter A = uj z /uj± is a measure of 
the aspect ratio between the axial confinement along the 
z direction and the radial confinement in the xy plane, 
with ui±_ = lu x — Lu y . The system is in the shape of a 
flat disk if A ^ 1 or a long cigar if A <C 1, including the 
limits A = 1 (spherically symmetric), A = (confined 
along the z axis), and A = oo (confined in the xy plane). 
The interaction between any two particles Vjk is taken 
to be a hard-sphere potential, modeled after the s-wave 
scattering lengths of the experimental systems, with an 
and <i22 being the diameters of the particles in species 1 
and 2, respectively, and a±2 — (an + 0,22) /2 being the 
potential range between two different particles. Masses 
mi and 777,2 and total numbers of particles iVi and N2 
are used for each of the two species, respectively, with 
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N = Ni + N2 being the total number of particles in the 
system. 

We performed simulations by carefully choosing the 
total number of time slices and the size of the Monte 
Carlo steps to ensure a good convergence of the physi- 
cal quantities calculated, including the density Pi(r) and 
the in-plane angular pair-correlation function T^(ip) — 
(8{lp — \ipi — <fk\)) f° r I ^ k. In order to examine the 
condensate profiles against the planar symmetry of the 
trap, we set the temperature much lower than the criti- 
cal temperature and projected all the particles into the 
xy plane when we take snapshots or calculate the aver- 
age radial density Pi{p) and angular correlation function 
Ti(ip), where p — y/ 'x 2 + y 2 is the in-plane radius and ip 
is the relative azimuthal angle. The density is normalized 

b y 2n Jo°° Pi(P)P d P = N i- 

In Fig.QJwe show three sets of simulation results for dif- 
ferent aspect ratio A under the same trapping frequency 
for both species at T/N 1 / 3 =0.1. It is clear that the 
symmetry of the Hamiltonian is preserved in the con- 
densate profiles when the trap is spherically symmetric, 
that is, A = 1. For both A = 4 and A = 16, the profiles 
are obviously asymmetric with a random orientation due 
to the spontaneous symmetry-breaking. The asymmetry 
becomes more severe and the correlation at small angle, 
as shown in Fig.^d), becomes much stronger as the ratio 
A increases. 

The symmetry-breaking in one species compensates 
the other because the center of mass of the system, un- 
der the given conditions, remains at the center of the 
trap, resulting in a stronger deformation in the lighter 
species. We performed one more simulation with A < 1 
and did not see any symmetry-breaking of the profiles in 
the xy plane. This is consistent with the findings of pre- 
vious mean-field calculations 0, 0] based on the Gross- 
Pitaevskii equation, which concluded that the symmetry 
is broken in the direction of the weakest trapping fre- 
quency. 

Let us now examine how the system responses to the 
change of the ratio 77 = aix/a-zi-, which is a measure of the 
imbalance of the particle sizes between the two species. 
In Fig. |2 we show three sets of snapshots with different ij 
while having other parameters mijm\ = 4, N\ — N2 — 
100, A = 16, and a n = 0.15 fixed. 

For the case of r\ = 1/3, the system shows a condensed 
core of species 1, surrounded by an outer ring of species 
2. Both condensates appear to have the planar symme- 
try of the trap. When rj is increased to 1/2, the system 
undergoes a quantum phase transition from two separate 
condensates to a binary mixture with each breaking the 
planar symmetry. This is evident from both snapshots 
and the correlation functions shown. The correlation 
at small angle increases when the symmetry is broken. 
Note that the condensate of species 1 is expected more 
off-centered ( 4 times) because the center of mass of the 
whole system must remain at the trap center. One can 
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FIG. 1: Snapshots of multi-exposures of the particles viewed 
along the z direction for each species in a double-condensate 
mixture, confined in a trap with 7712/mi = 4, Ni = N2 — 100, 
an = 0.2, and 022 = 0.4, for different aspect ratio A = u} z /uj±: 
(a) A = 1; (b) A = 4; and (c) A = 16. The corresponding in- 
plane angular pair-correlation functions Fi(tp) are shown in 
(d). The planar symmetry in each of the two condensates is 
broken if A 7^ 1; the asymmetry becomes more severe when A 
is further away from 1, especially for the lighter species. 



see from the snapshots that the binary phase is formed 
from the outward movement of species 1 plus the inward 
movement of species 2. Therefore, increasing r/ further 
can result in exchanging the roles of the two species. This 
is precisely what happens when rj is increased to 1 — the 
system shows a condensed core of species 2, surrounded 
by an outer ring of species 1. The correlation at small 
angle decreases and the planar symmetry of the Hamilto- 
nian gets restored in the condensate profiles in the same 
time. 

To investigate the temperature dependence of asym- 
metry of each condensate, we divide the trap into two 
equal halves by a plane through the z axis, with the low- 
density side containing a minimum number of particles 
(-/Vl) and the high-density side containing a maximum 
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FIG. 2: Snapshots of multi-exposures of the particles viewed 
along the z direction for each species in a double-condensate 
mixture, confined in a trap with 7712/mi = 4, Ni = N% = 100, 
A = uj z /lo± — 16, and an = 0.15 fixed, for different ratio 
77 = 011/022: (a) 77 = 1/3; (b) 77 = 1/2; and (c) 77 = 1. 
The corresponding in-plane angular pair-correlation functions 
are shown in (d). The planar symmetry of the particle 
distributions is broken first and then restored when the ratio 
is varied monotonically. 



number of particles (-/Vh). Then we examine the ratio 
£ = TVl/^Vh, which is shown in Fig. [3] for a trap with 
77 = 1/2, m 2 /m 1 = 4, ATi = N 2 = 100, and A = 16 fixed. 

A general feature emerges from all these systems: The 
asymmetry is depleted in both the extremely low tem- 
perature region and high temperature region. The causes 
are however different. At the extremely low temperature, 
the strong coherence in each of the condensates forces 
the particles to go to the lowest single-particle states, 
and hence to form a nearly symmetric distribution. This 
is extremely evident from the systems with a weaker 
interaction, for example, in cases with an = 0.1 and 
an = 0.15, where the symmetry is nearly completely re- 
stored when temperature is lower than 0.057V 1 / 3 . At the 
high temperature, particles gain more thermal (kinetic) 
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FIG. 3: Temperature dependence of the condensate profile 
asymmetries, measured from the ratio £ = N-l/Nh, where 
TVl is the number of the particles in the low-density half and 
TVh is the number of particles in the high-density half, with 
77 = 011/022 = 1/2, m 2 /m 1 = 4, Ni = N 2 = 100, and A = 16 
fixed. The most severe asymmetry for each case happens at a 
temperature that decreases with the interaction (hard-sphere 
diameter) . 



energy and can therefore spread out more in order to 
distribute more symmetrically. More detailed analysis of 
the data also show that for weaker interaction cases, the 
particles from one condensate penetrate into the other 
condensate when temperature is increased, whereas for 
the stronger interaction cases, the particles from one con- 
densate go around to form an outer ring. 

The optimal temperature, at which the system shows 
a minimum Nl/N-r (maximum asymmetry) in Fig. [31 
decreases with the interaction range. This is consistent 
with our earlier findings 17] that the many-body corre- 
lation effect, originated from the strong interaction, is in 
competition with the coherence and it is more difficult to 
restore the symmetry in the systems with a strong inter- 
action. Another distinctive feature obtained here is that 
the minimum ratio of Nl/Nh is smaller if the interac- 
tion range is greater. In other words, the system with a 
stronger interaction can reach a higher asymmetry. 

In addition to the interaction and temperature, the 
structure of a double-condensate mixture also depends 
on the total number of particles in the system. In Fig.^ 
we show the results of three different N 2 while keeping 
Ni = 100, m 2 /m 1 = 4, an = 0.2, a 22 = 0.4, and A = 16 
fixed. 

From the correlation functions shown for N 2 = Ni = 
100, the particles from species 1 occupy one side of the 
trap with a significant asymmetry while the particles 
from species 2 stay at the center of trap with a smaller 
asymmetry, also seen in Fig. 12(c). Quantitatively, one 
expects that the center of mass of species 1 is about four 
times away from the trap center of that of species 2 in 
this case because the center of mass of the system is at 
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FIG. 4: (a) Density profiles Pi(p) and (b) angular pair- 
correlation functions Ti (ip) of the particles confined in a trap 
with iVi = 100, m 2 /mi = 4, an = 0.2, a 2 2 = 0.4, and 
A = 16 fixed, for N2 = 100 (circles), 200 (triangles), and 
300 (crosses). As N2 increases, the lighter particles, while 
maintaining certain asymmetry, contract and are eventually 
surrounded by the heavier particles that always appear to be 
distributed symmetrically and are pushed out of the center. 

the center of the trap. When N2 is increased to 200, 
the particles in species 2 squeeze each other outward but 
push some of the particles in species 1 into the central 
area and others further off the center. Finally when N2 is 
increased to 300, particles in species 1 all move to the cen- 
tral area; both condensates recover the planar symmetry 
of the trap with particles in species 2 forming an outer 
ring, surrounding completely the particles from species 
1. However when TVi is increased from 100 to 300 with 
N2 kept to be 100, the asymmetries return with each 
condensate occupying one side of the trap. 

So far we have shown the results with both the species 
under the same trapping frequency, which means a dif- 
ferent external potential for each species if m 2 /mi ^ 1. 
In order to examine the influence of external potential, 
we can either choose m^/nii = 1 with the same trap- 
ping frequency or set 7711/7712 = uj\jiJ\. Such simulations 
performed for those systems depicted in Fig. show no 
broken symmetry as the aspect ratio A is increased from 
1 to 16, although phase separation still occurs along p di- 
rection. This is similar to the case of our previous work 
on the double-condensate mixture in a spherically sym- 
metric trap we a l so nn( i that the particles with a 
larger interaction range are in the outer ring while the 
other condensate stays at the central area. 

The external potential plays an important role in the 
development of asymmetry in the systems studied. For 
a system with 1712/ mi = 4 under the same trapping 
frequency, the condensates can undergo a symmetry- 
to-asymmetry transition and then an asymmetry-to- 
symmetry transition, as shown in Fig. [2J However, for 
a system with 1112/ mi = 1, still under the same trapping 
frequency, no significant asymmetry can be developed. 
The system with 011/022 = 0.5 can have a stable, phase- 
separated configuration similar to that of Fig. [21 a), which 
is also a stable configuration obtained in a mean-field 



approach || with ai 2 /(aiia 2 2) 1 ^ 2 ~ 1-06. Furthermore, 
when an/a22 is increased to 1, the two separated con- 
densates penetrate each other and form a binary mixture 
with a slight asymmetry — the double condensates are 
mirror images of each other because of the interchange- 
able nature of the two species under the given conditions. 

In summary, our simulations show that asymmetry can 
occur in a double-condensate mixture confined in a disk- 
shaped trap. The asymmetry has a strong dependence 
on the aspect ratio of the axial confinement along the 
z direction and the radial confinement in the xy plane, 
the ratio of the interaction ranges 011/022, and the tem- 
perature. Furthermore, the total numbers of particles in 
the two species and significant difference in the external 
potentials can also affect the structures of the conden- 
sates. The novel behavior of a double-condensate mix- 
ture elucidated in this work can be verified experimen- 
tally with the setup that has created the mixtures and it 
is extremely interesting to see such phenomena involving 
quantum phase transition and reentry of symmetry of a 
model Hamiltonian realized in a laboratory. 

This work is supported in part by DOE under Co- 
operative Agreement DE-FC52-01NV14049 and by NSF 
under Cooperative Agreement ACI-9619020 through the 
computing resources provided by the National Partner- 
ship for Advanced Computational Infrastructure at the 
University of Michigan Center for Advanced Computing. 



[1] C.J. Myatt, E.A. Burt, R.W. Christ, E.A. Cornell, and 
C.E. Wieman, Phys. Rev. Lett. 78, 586 (1997); D.S. Hall, 
M.R. Matthews, J.R. Ensher, C.E. Wieman, and E.A. 
Cornell, ibid. 81, 1539 (1998). 

[2] C.A. Regal, C. Ticknor, J.L. Bohn, and D.S. Jin, Nature 
424, 47 (2003). 

[3] J. Herbig, T. Kraemer, M. Mark, T. Weber, C. Chin, 
H.-C. Nagerl, R. Grimm, Science 301, 1510 (2003). 

[4] HP. Buchler and G. Blatter, Phys. Rev. Lett. 91, 130404 
(2003). 

[5] R.A. Duine and H.T.C. Stoof, Phys. Rev. Lett. 91, 
150405 (2003). 

[6] S. Gupta, Z. Hadzibabic, M.W. Zwierlein, C.A. Stan, 
K. Diechmann, C.H. Schunck, E.G.M. van Kempen, B.J. 
Verhaar, and W. Ketterle, Science 300, 1723 (2003). 

[7] B.D. Esry and C.H. Greene, Phys. Rev. A 59, 1457 
(1999). 

[8] A. A. Svidzinsky and S.T. Chui, Phys. Rev. A 67, 053608 
(2003). 

[9] P. Ao and S.T. Chui, Phys. Rev. A 58, 4836 (1998). 
[10] P. Ohberg and S. Stenholm, Phys. Rev. A 57, 1272 
(1998). 

[11] P. Ohberg, Phys. Rev. A 59, 634 (1999). 
[12] S.T. Chui and P. Ao, Phys. Rev. A 59, 1473 (1999). 
[13] J.G. Kim and E.K. Lee, Phys. Rev. E 65, 066201 (2002). 
[14] V.S. Shchesnovich, A.M. Kamchatnov, and R.A. 

Kraenkel, Phys. Rev. A 69, 033601 (2004). 
[15] T.L. Ho and V.B. Shenoy, Phys. Rev. Lett. 77, 3276 



5 



(1996). [18] S. Pearson, T. Pang, and C. Chen, Phys. Rev. A 58, 

[16] F. Riboli and M. Modugno, Phys. Rev. A 65, 063614 1485; 4796; and 4811 (1998). 

(2002). 

[17] H. Ma and T. Pang, cond-mat/0309004 



